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Abstract 

Uniformity of the probability measure of phase space is considered in the 
framework of classical equilibrium thermodynamics. For the canonical and 
the grand canonical ensembles, relations are given between the phase space 
uniformities and thermodynamic potentials, their fluctuations and correlations. 
For the binary system in the vicinity of the critical point the uniformity is in- 
terpreted in terms of temperature dependent rates of phases of well defined 
uniformities. Examples of a liquid-gas system and the mass spectrum of nu- 
clear fragments are presented. 

PACS: 05.70.-a, 64.60.-i 
Introduction 

Generalized entropies were proposed by Renyi [[p to characterize probability 
measures. Given a probability measure dp and a partitioning of a phase space into 
M(A) cells, {Ai}^^', each of the same volume A, the Renyi entropies I q (q 6 R 1 ) 
are defined 

'AlnLS A y, q*1 

Iq = < (1) 



V M(A) , i 

2_i=i Pi ln Pi. q = 1 



where 



dp (2) 

and sums run over cells for which pi ^ 0. It has been found later by many authors 
[j|] that for investigation of singular and multifractal measures a useful quantity is a 
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generalized dimension D q , related to the Renyi entropy I q 



D q =-lim-\. (3) 

a^oo In A 

The Renyi dimension D q is always positive and decreases with q, so that there al- 
ways exists Dqo = liiriq^oo D q . For positive q, the D q and I q is mostly sensitive to 
these regions of the phase space where the measure is concentrated and for negative 
q - where it is rarified. 

In order to find a universal description of probabilty distributions for multifrac- 
tals, another q -dependent quantity, called uniformity, was constructed by Beck in 
ref. [|3|] from Renyi dimensions or entropies 

_ D q + q(q-1)D q 
7q " D q + (q-1)D q W 

where ' stands for derivative over q. For all q, < y q < 1 . It can be found solving 
eq. (Q) that for the most uniform probabily density y q = 1 and 

D q = D 0C (5) 

does not depend on q. The minimum uniformity y q = corresponds to 

1 

1 -"l/q 

The derivative of the uniformity 



D q = i tttDoo- (6) 



, ^ D q [2D q -(1-q)D^] 
y ^ (q " 1] [D q -(1-q)D q P (7) 

vanish for either q = 1 , if it is differentiable there, or for 

D q = D 00 + ^ (8) 

1 - q 

The uniformity was introduced in order to characterize random variables on 
events generated by formal systems. It seems interesting to apply this concept to 
random variables related to physical systems in various conditions and represented 
by variables of thermodynamical ensembles. Since the probability measures for 
such ensambles are well defined and they determine both Renyi dimensions and all 
traditional thermodynamic potentials, the relations between these quantities should 
be easy to find. On this way one could expect to gain some understanding on how 
the structure of the phase space is reflected in such global characteristics of the sys- 
tem as its thermodynamic functions. Moreover, because of sensitivity of D q and 
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Yq, controlled by the parameter q, to singularities of the measure, their relation to 
fluctuations of some measurable quantities, as internal energy or numbers of parti- 
cles, or correlations of those, has to be straightforward. It can be also expected that 
two control parameters, the q for Renyi entropies and the temperature for thermo- 
dynamic potentials, should to some extent play similar roles in the formalism, the 
latter determining the strength of thermal fluctuations and the former steering the 
sensitivity of D q and Yq to them. 



The Canonical Ensemble 

Consider the canonical ensemble with the probability measure 



where E stands for energy and 



Z(|3) 



dn (10) 



is the partition function, where the integral is taken over the whole phase space, 
[3 = 1 /T is inverse temperature and the Boltzmann constant is assumed equal one. 
Using (^) we find for A — > 



F(T/q)-F(T) 



q^i 



T(1/q-1) ' 

D q lnA = { (11) 
_™ = -S(T), q = l 

where 

F(T) = -TlnZ(T) (12) 

is free energy. Relation similar to ([TT]) , between Renyi entropies and free energy, 
was found by Tel in different context in ref. [Q]. It is seen from eq. ([[]]) that q is 
temperature scaling parameter. For q ^ 1 the Renyi entropy I q is equal to the finite 
difference ratio of F and T and for q — > 1 it becomes the usual relation between F 
and the thermodynamic entropy S. Using (JTT|), the uniformity (|]) for the canonical 
ensemble is found 

J_ 1 F(T)-F(T/q) 

Y q q TS(T/q) ' 1 } 

The maximum uniformity Yq = 1 corresponds to 

qF(T/q) = F(T) = — TS max (14) 
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where S max is the maximum entropy. For the minimum uniformity, y q = 0, the free 
energy F does not depend on q and 



F(T/q) = F(T) = -TS min , q>1 (15) 

where S m ; n is the minimum entropy. 

Calculating the first and the second derivatives of D q and substituting to eq. (|7|) 
one finds 

^'^WFW™' 11 - (16) 

where V(E) is the variance of the energy E and U = (E) is internal energy, or 
rewriting it in terms of the specific heat Cy(T/q) = V[E(T/q)]/(T/q) 2 , using eq. 
([□) and formula F = U - TS 

3, q > 1 

y'aVn = { (17) 

Cy(T/q)(1-q)S(T) i- 
[{l-q)S(T)-S(T/q)F» H ~> 1 



The Grand Canonical Ensemble 

For the grand canonical ensemble the probability measure depends in addition 
on the number of particles N 

d^ N = —J— e p[ ^ N - EN] dE N (18) 

where \i c is the chemical potential and E is the grand canonical sum 



N=0 



du- N . (19) 



The integral is taken over the whole N -particle phase space. 
The Renyi dimensions in this case are (A — > 0) 

( (N(T/q))-(N(T)) , , 

M-c T(1/q-1) ' ^ T 1 

D q lnA = I (20) 
[ -^m = _s(T), q = l 

where the average number of particles is given by 

(N(T))=T ain|a if d 

o\i c 
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and the uniformity is equal to 

1 _ 1 (N(T))-(N(T/q)) 

y q "q + ^ C TS(T/q) " (22) 

Analogously to the canonical ensamble, the condition for the maximum uniformity 
is 

q(N(T/q)) = (N(T)) = -S max (23) 

and for the minimum uniformity is 

(N(T/q)} = (N(T)) = — S^, q>1. (24) 

M-c 

Calculating the first and the second derivatives of D q one finds the derivative of 

Vm = 1 (N(T/q))-(N(T)> 

yql ' T/q [U(T/q) - M(N(T/q)) - (N(T)»P 

{V[M- c (N(T/q))]+V[E(T/q)]-2^ c cov[(N(T/q)),E(T/q)]}. (25) 

In case of constant volume and using eq. flZU|), the formula can be rewritten in the 
vicinity of q = 1 in terms of the specific heat Cy and isothermal compressibility 
K T (T/q)=V[(N(T/q)>]/p(N(T/q)> 

f 0, q>1 



y'm = 



(l-q/T)S(T) (f)f ^ 
IT/qS(T/q)+Mc<N(T)>F X 

{p c p KT (T/q)(N(T/q)) + (T/q) 2 c v (T/q)- 

l 2p c cov[(N(T/q)),E(T/q)]} ) q -> 1" 



Applications and Discussion 

The uniformity parameter y q was originally proposed [^]] to measure distance 
from the two extremes: the chaotic phase (y q = 1) and condensed phase (y q = 0). 
We show in the following that for the system consisting of subsystems in thermal 
equilibrium, the overall uniformity is the weighted average of uniformities of sub- 
systems with weights given by particle numbers ratios for subsystems. 

Consider the two-phase system with average numbers of particles (Ni(2)) and 
the uniformities of pure phases Yi[2> I n order to find the condition for Yi[2} to De 
independent of q one solves eq. (Q) 

Pi(2) q + q(q-1)D; (2]q 

L>1(2)q + lq - I J U l(2)q 
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with respect to D , 



D 



K2)q 



V' 



K2)q 



D l(2)oo 
D l(2)oo 



g -7K2) 

q-1 

Tl(2) ~ 1 
(q-1) 2 



(28) 



Using eq. ( |20| ) and the fact that chemical potentials in subsystems are equal in 
equilibrium, one finds 



D q = D lq + D 



2q 



and, using ( |28| ) 



where 



Tq(T, \i c ) = Ci(T, M- C )y 1 + c 2 (T, [i c )y 2 , 



C 1{2) 



(Nl(2)} 



(N 1 ) + (N 2 ) 



(29) 



(30) 



(3D 



are concentrations of particles in subsystems 1 and 2. Generalization to more than 
two subsystems is straightforward. 

Example 1: classical liquid-gas system 

Consider any classical system of two coexisting phases in thermal equilibrium 
near the critical point. In the vicinity of critical temperature T c , phase densities 
follow the power-law dependence on temperature [R 



P2 - Pi = const • (T c - T) 



T< T r 



(32) 



where <x is the critical exponent (e.g. for the liquid-gas system its value is between 
0.33 and 0.36). Assuming the average total number of particles (Ni + N 2 ) inde- 
pendent of temperature and the volume of the system to be constant, from eq. ( |32| ) 
follows that 



c, = a + b(i-I 

C 2 = 1 - Ci 



(33) 



where A and B are constants. As discussed before, the parameter q plays the role 
of the temperature scaling parameter. Thus the temperature evolution in eq. ( |33| ) 
can be rewritten in terms of q = T/T c Q. Then q = 1 corresponds to the critical 

'The choice of T c as the reference is arbitrary but convenient. For another temperature scale one 
has to change the boundary conditions in eqns. ( [34] ) 
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temperature and domains q > 1 and < q < 1 correspond to pure uniform (gas) 
and mixed phases, whereas the limit q = should represent pure condensed (liquid) 
phase. For temperature T much lower than T c there is only condensed phase in the 
system and a weak dependence of density on T can be neglected. Therefore, for 
q = T/T < 1 , one finds the values of A and B from the limits 

lim Yq = yt 

q->1- 

lim Yq = 72- (34) 

q-> qo 



Hence, from (33) 



1 /1-q\ a 



1 - Y2/T1 V 1 - qo 
J f 1-q 

i -72/71 v 1 - qo 



Cl 



c 2 = , -) 05) 



for q < q < 1 and 



Cl(2) = 

c 2( i] = 1 (36) 

for q < qo(> 1 )• The Yq does not have to be differentiable at q = 1 . 

The uniformity Yq for <x = 0.34 and Ci 2 given by ( |3"5| ) is shown in fig. 1. Since 
the power-law dependence of densities difference on temperature is valid for the 
near-to-critical region, this picture is realistic only for T — > T c or q — > 1 (see broken 
horizontal axis). However, if temperature dependence of phase densities is known 
for any temperature, the same line of reasoning can be followed to find the overall 
uniformity for the mixture of phases in wide range of temperatures. For some mod- 
els, e.g. the power-law probability density mentioned in ref. [Q], the uniformities 
of pure phases can be exactly calculated. In other cases, provided the probability 
measure is known, these pure uniformities can be also determined by direct, al- 
though sometimes complex and time consuming, computation. On the other hand, 
if the phase space probability density in the mixed-phase domain is measured ex- 
perimentally and the uniformities of pure phases are known from elsewhere, or can 
be calculated, this theoretical scheme can be used to determine the ratio of phases 
at given temperature. 

Example 2: discrete power-law energy spectrum 

Let us discuss the system of discrete energy spectrum, where the energy E of 
the system is a random variable taking positive integer values and exhibiting the 
power-law probability measure 

d\i E = -— 6(E-m)dE, m = 1,2,... (37) 

Cm 
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1 



q=T/i c 



Figure 1: The uniformity y q for the two-phase system, assuming the near-to-critical 
power-law behaviour of densities difference as functions of temperature. Critical 
exponent <x = 0.34 is assumed. 

where C(t) = Y-°^\ i"n-~ T , (t > 1 ) is Riemann zeta function. Expected value of 
the energy is then 



and exists only for t > 2. Such statistical ensemble would be interesting per se and 
similar study of its properties in terms of Renyi entropy and Beck uniformity, as for 
the canonical and grand canonical ensembles, could be performed. 

In order to stay within the framework of equilibrium thermodynamics and use 
the concept of equilibrium temperature, consider the canonical case. The statistical 
sum is equal to 




C(T-l) 



(38) 



oo 



Z(T) = ^e- m/J 



m=l 



1 



(39) 



Internal energy is equal to 



U(T) = 



91nZ(T) 
3(VT) 



1 - e-VT 



(40) 
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and by requiring it to be equal to (E) we find the equiliblium temperature 

1 

ln(1-~1/ i 
and the equilibrium probability measure 



T eq = (41) 



d^ E = 1 6(E-m)dE, m = 1,2,.... (42) 

For this probability density we find (A — > 0) 

-^ln[<E)<-«E)-1)«], q^l 
D q (T eq )lnA = <( (43) 

(1-(E))ln((E)-l) + (E)ln(M), q = l 



and 



[{E)q-((E)-1)q]ln[(E)q-({E)-1) q ] n /I 

4 (E)iln(E)-((E>-1)iln({t)-1) ' H 7= 1 

Tq(T e c,) = <; (44) 

1, q = l 



Table 1: Power-Zaw exponents, temperatures, mean energies and uniformities for 
nuclear fragmentation data. 



T 


(£> 


T(MeV) 


q = T/T c 




4.1 


1.10 


6.0 


0.50 


0.29 ±0.05 


3.8 


1.14 


6.2 


0.52 


0.36 ±0.07 


3.7 


1.15 


6.7 


0.56 


0.49 ± 0.08 


3.0 


1.37 


7.2 


0.60 


0.65 ±0.11 


2.8 


1.51 


7.7 


0.64 


0.75 ±0.12 


3.2 


1.28 


8.2 


0.68 


0.79 ±0.12 


2.6 


1.74 


8.3 


0.69 


0.83 ±0.13 


2.4 


2.21 


14.0 


1.16 


0.98 ±0.12 


2.9 


1.43 


14.5 


1.20 


0.96 ±0.11 


2.9 


1.43 


15.0 


1.25 


0.95 ±0.11 



As an example of a real system exhibiting the power-law energy spectrum we 
employ the nuclear mass spectrum of heavy nuclei fragmentation. Taking simplify- 
ing assumptions that the whole energy is identical to the fragment mass, neglecting 
the surface energy, the binding energy and finite size effects and assuming that all 
states are distinguishable, we use the power-law exponents determined experimen- 
taly for low mass fragments and refered to in ref. The power-law behaviour is 
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Figure 2: The uniformity y q for experimental data from nuclear fragmentation. The 
dashed line represents analytical uniformity for x = 4.1 and the dotted one for 

T = 2.1. 



also characteristic for liquid cluster size distribution near the liquid-gas transition 
point in real systems and in percolation models. It has to be noted however, that T 
exponent obtained for nuclear systems disassembly [Q| is not a universal critical ex- 
ponent but an apparent exponent absorbing some temperature dependence, usually 
unknown. Nuclear data were analysed in ref. [§]] in the framework of the conden- 
sation theory and the critical temperature T c = 12.0 ± 0.2 MeV was found. For 
each data sample exhibiting the power-law mass distribution the temperature was 
determined from either the slope of the energy spectrum at 90° or from the Fermi 
gas model or from the moving source model. We identify this temperature with the 
equilibrium temperature ([41]), define q = T/T c and calculate y q (T) from eq. (fHj). 
The numbers are given in tab. 1 and the uniformity as a function of q in fig. 2. For 
the calculation of errors of y q the values for errors of T c , T and T are taken from ref. 
[§]. The curves in fig. 2 represent uniformities calculated using analytical formula 
(H4j) for t = 2.4 and t = 4.1 . As all experimental points stay within the belt deter- 
mined by these curves we conclude that the temperature evolution given by ( fPD is 
consistent with data. The data show that uniformities increase with temperature in 
the mixed-phase region and are consistent with 1 above critical temperature. Below 
T = 4 MeV the contribution of the non-uniform phase is below 1 0%. 
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